Orbit and Optics Improvement by Evaluating the Nonlinear BPM Response in CESR 
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We present an improved system for orbit and betatron phase measurement utiUzing nonUnear 
models of BPM pickup response. We first describe the calculation of the BPM pickup signals as 
nonlinear functions of beam position using Green's reciprocity theorem with a two-dimensional for- 
malism. We then describe the incorporation of these calculations in our beam position measurements 
by inverting the nonlinear functions, giving us beam position as a function of the pickup signals, and 
how this is also used to improve our measurement of the betatron phase advance. Measurements 
are presented comparing this system with the linearized pickup response used historically at CESR. 



I. INTRODUCTION 
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CESR measures beam position and betatron phase 
with approximately one hundred beam position moni- 
tors (BPMs) distributed around the storage ring. Each 
BPM consists of four button-type electrodes mounted 
flush with, and electrically isolated from, the surface of 
the beam pipe. A moving particle bunch induces charge 
on the beam pipe walls and on the surface of each button, 
which one can describe as image currents or as surface 
charge due to the transverse component of the bunch's 
electric field [l|. 

The BPM buttons are connected to electronics that 
process and record signals which are a function of the dis- 
tance between the button and the passing bunch. The 
four signals from each BPM are used to determine the 
beam position and betatron phase advance. At many 
accelerators, the button signals' nonlinear dependence 
on the beam position is linearized for simplicity. Before 
the improvements described here, this approach was also 
used in CESR. 

Our efforts to improve the beam position measure- 
ments by including the nonlinear BPM response is mo- 
tivated by CESR's pretzel orbits, where electron and 
positron beams avoid parasitic collisions by following sep- 
arate paths with large displacements from the central 
axis of the beam pipe. The linearized methods are not 
reliable for such large amplitudes, and have made accu- 
rate beam position and betatron phase measurements at 
CESR impossible under colliding beam conditions. We 
will illustrate those shortcomings and present measure- 
ments demonstrating improvement by using the nonlin- 
ear models. 



II. BACKGROUND 



Many accelerators, including CESR, have traditionally 
assumed a linear relationship between the beam posi- 
tion and the BPM button signals. Given four signals 
Si {i = I, ... ,4:) from buttons arranged as in Fig. ^ the 



5-2 



FIG. 1: Arrangement of buttons in CESR arc BPMs 



transverse beam position is given approximately by 
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where k^^y are scale factors set by the geometry of each 
BPM type. This evaluation of BPM signals is often called 
the difference- over- sum method. Equations (|1I2|I pro- 
vide an estimation of the bunch position in relatively few 
arithmetic operations. 

Because analytical approaches to determining the fac- 
tors kx.y make drastic approximations to the BPM geom- 
etry, we have tried to measure the factors experimentally 
at CESR through a variety of techniques summarized in 
Table m Those include translating a section of the beam 
pipe containing the BPM with precision actuators, us- 
ing a test stand with a movable antenna to simulate the 
beam, and using the known value of the dispersion while 
changing the beam energy in dispersive regions |3|. 

But precise knowledge of kx,y is of limited benefit, since 
Eqs. (|1I2|I yield only the linear part of the signal depen- 
dence for bunches near the center of the BPM. In the next 
section, we describe our technique for accurately calcu- 
lating button signals, but let us first use those results to 
illustrate the limitation of the linearized formulae. 

In Fig. 12 the four button signals were calculated nu- 
merically for a beam at different horizontal displace- 
ments. The four signals were combined according to 
Eq. ^ and plotted against the known horizontal dis- 
placement. The slope near the origin gives k~^, but the 



Method 


kx (mm) 


ky (mm) 


20 MHz antenna 


25.58±.33 


20.58±.43 


Dispersion (1990) 


26.3 




Dispersion (1991) 


27.4±.6 




Moving beam pipe (e^) 


26.82±.25 


19.96±.ll 


Moving beam pipe (e~) 


27.14±.54 


20.48±.19 


2D Poisson model 


26.2 


19.6 



TABLE I: Measured scale factors for CESR arc BPMs. 
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FIG. 2: Linearized horizontal position measurement in an arc 
BPM. At a typical pretzel amplitude of 1.5 cm, the linearized 
formula shows significant disagreement with a linear fit. 



linear relationship breaks down noticeably at approxi- 
mately 1 cm, and beyond 2 cm the relationship fails com- 
pletely. Because pretzel orbits in CESR are typically as 
large as 1.5 cm, this is precisely what has hindered accu- 
rate measurements under colliding beam conditions until 
the improvements described in this paper were imple- 
mented. 
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FIG. 3: Linearized map distortion in CESR interaction-point 
BPM with approximately circular cross-section. 
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FIG. 4: Linearized map distortion in CESR arc BPM with 
approximately elliptical cross-section. 



The problem of this nonlinearity is made even more ev- 
ident in two dimensions. Figures I3I4I show a regular grid 
of (x, y) points and the mapping of those same points un- 
der Eq. (|1I2|I . Both BPMs show the characteristic pin- 
cushion distortion, which increases with distance from 
the origin. 

In CESR, betatron phase measurements also rely on 
a related assumption about the linearity of the button 
signals. The betatron phase is measured by shaking the 
beam at a sideband of the betatron frequency. For each 
detector, the phase for each button is calculated by elec- 
tronically comparing AC signal on that button to the 
phase of the shaking. From the individual horizontal 
(vertical) button phases O^a i&v,i), the horizontal and ver- 
tical betatron phase is calculated by 






r,^tll,h 



93,H (3) 

'^-, (4) 

where A^^v are real constants that are not used fur- 
ther 13 . Other than the minus signs which account for 
the assumption that the beam is shaking between the 
pairs of buttons, this is simply an averaging of button 
phases represented as complex vectors. 

When the horizontal orbit amplitude is large, the beam 
begins to shake underneath the buttons, and the rela- 
tionship between the beam motion and the button signal 
becomes complicated. In such cases, some of the buttons 
may report an inaccurate phase, and averaging them with 
the rest corrupts the final answer. We will show how our 
nonlinear models can improve not only beam position 
measurements, but these measurements as well. 



III. AN IMPROVED SYSTEM FOR POSITION 
AND PHASE MEASUREMENT 

In order to overcome the limitations described, a new 
system has been implemented with two major compo- 
nents: realistic numerical models of the button response, 



and an efficient algorithm for inverting the model to yield 
beam position. 



A. Numerical Calculation of BPM Response 

For accurate beam position measurements, a function 
is required that expresses the bunch location (x, y) as a 
nonlinear function of the button signals. Since the four 
button signals lead to two coordinates (and a scale fac- 
tor), the problem is over-constrained, and this function 
cannot be obtained directly. The inverse (button signals 
from beam position), however, is readily obtainable by 
standard numerical techniques. 

The most accurate and direct method of numerical so- 
lution is to simulate the bunch in a three-dimensional 
BPM, calculating the electromagnetic fields, and from 
them, the charge on the buttons. The simulation could 
be repeated for different beam locations, and the fields 
recalculated until enough solutions were accumulated to 
describe the behavior over the entire BPM. However, this 
is very computationally intensive, and can be avoided by 
the methods that follow. 



1. Two-Dimensional Approximation 

For ultra-relativistic bunches in a beam pipe with 
constant cross section, the approximate electromagnetic 
fields can be calculated using a two dimensional formal- 
ism IJ, Q ■ Assuming the bunch has negligible transverse 
extent, the charge distribution of the bunch may be writ- 
ten, in the lab frame, as 



= 5{v - To) ^ pk cos {k(z - vt)) 



(5) 



where the longitudinal dependence in z has been written 
as a Fourier expansion. Transforming to the reference 
frame of the bunch, the charge density and electric po- 
tential are written 
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cos(fcz*/7). 



(6) 

(7) 



We write Poisson's equation V^^* — p* in the bimch 
frame as 



Vi - — ) $(r)0fe - 5{v - vo)pk (8) 



where V^ is the two-dimensional transverse Laplacian. 
For bunches with length ai without appreciable longitu- 
dinal substructure, pk is only relevant for fc < — . The 
characteristic distance over which <i>(r) changes is the 
diameter a of the beam-pipe so that the order of mag- 
nitude estimate |V^<i>(r)| « -^\^\ can be made. For 



sufficiently long bunches and sufficiently large values of 
7, the relevant values of k/^ can be neglected, i.e. when 
72 "^ (^)^ ^'^'^ ^^^ solution is described by the two di- 



mensional, electrostatic case 



Vl$(r) = ^<5(r-ro) 

(Pk 



(9) 



Since we only need $(r) up to a multiplicative factor, we 
don't worry about the constant coefficients on the right- 
hand side. 



2. Green's Reciprocity Theorem 

Rather than perform a separate calculation of the but- 
ton signals for many beam positions, we use Green's reci- 
procity theorem to calculate the button signals for all 
(x, y) inside the BPM with a single numerical calcula- 
tion. This theorem states that the surface charge a on 
a button due to a test charge at {x, y) is proportional to 
the potential at that same position when the test charge 
is absent and the button is excited by a potential V. 

Suppose we have two scalar functions (pi and (/)2 in a 
volume V bounded by a surface S. We form the vector 
field 

A = (/)iV(/.2 (10) 

for which the divergence theorem guarantees 
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Manipulating the integrands gives 

V-(0iV02) = (V0i)-(V02) + 0lV202 (12) 

A-fi = ^iV02-ft-0i^ (13) 

on 

where n is a unit vector normal to the surface and point- 
ing out of the volume of integration, and d /dn indicates 
differentiation with respect that direction. Equation Ijllf) 
yields 

[(V0i) • (V02) + 0iV202] dV = j 4>i^da. (14) 

If we interchange 0i and 02 and subtract the result from 
Eq. 114|) , we can eliminate the first term in the integrand 
of the left hand side. This gives 



[01 V^ 



i]dV= j> 



dn 



dn 



da. 

(15) 

Taking the (j)i to be potentials for volume charge den- 
sity Pi and surface charge density ai leads to Green's 
reciprocity theorem: 



1P2 dV 



11CT2 da ~ (f>2Pi dV 
s Jv 



2cri da 



(16) 



where we have used V^0 = —p and d(j)/dn = a (recall 
that n points into the conducting surface). 

Connecting this result to the case of a BPM, imag- 
ine 01 corresponds to the potential when a single button 
is excited with a potential V and all other surfaces are 
grounded. We can calculate the potential 4>i{x, y) by nu- 
merical solution of Laplace's equation. For the second 
potential (/)2, we ground all surfaces and put a charge 
distribution p2{x,y) inside the BPM. 

We plug the two cases into Eq. (|16|l and observe that 
the third integral vanishes because there is no volume 
charge for the first case (pi = in y ) . The fourth integral 
vanishes because we grounded the beam pipe and the 
buttons (02 = on 5) . Since V can be pulled out of the 
second integral, what remains is just the total charge on 
the button, labeled qi,, giving 



<Pi{x,y)p2{x,y)dV = -Vqt. 



(17) 



If p2 is a point charge q located at {xo,yo)i then the 
integral in Eq. I|17|l picks out the value (pi{xQ,yQ). We 
arrive at the final relation 
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q(f'ixo,yo) 
V 



(18) 



remembering that (j){x, y) and V refer to the two different 
configurations. 

Therefore, since the signal on a button is proportional 
to the induced surface charge on that button g^, (t>{xo,yo) 
is the solution to the problem of calculating the button 
signal, up to a multiplicative constant, as a function of 
the bunch location. 



3. Calculation of the Button Signals 

Based on the previous arguments, we use POISSON 
to solve the boundary value problem for 0(a;, y). For the 
two-dimensional boundary, we take a slice at the longitu- 
dinal midplane of each BPM. The first button is set (ar- 
bitrarily) at 10 volts and all other surfaces are grounded. 
POISSON generates a mesh inside the boundary, com- 
putes the solution to Laplace's equation on the mesh, 
and stores the result at regular grid points in an output 
file. 

CESR BPMs have multiple geometric symmetries, so 
the signals (j)i{x,y) on the other three buttons are just 
reflections or rotations of the coordinates for the excited 
button in the first calculation of 0i(x,y). To compute 
(j)i{x,y) between grid points, we use bicubic interpolat- 
ing polynomials, which are stored for quick subsequent 
evaluation. 



B. Realtime Inversion 

For beam position measurements, we start with button 
signals Si and seek the location (a;, y) of the beam. The 



result 4'i{x, y) from the POiSSOiV calculation must be in- 
verted, and since we have four constraints (four buttons) 
and three parameters (position {x,y), and a scale factor) 
we proceed by fitting the calculated button signals to the 
measured signals. We minimize the merit function 






{qMx,y) ~ Sif 



crt 



(19) 



where (j)i{x,y) is the signal on the i"^ button and the ai 
are the uncertainties in the measured signals (which we 
take to be the same for all four buttons). The factor q is 
proportional to the beam current and could be used for 
beam loss studies. 

Minimization is performed via the Levenberg- 
Marquardt method provided in Numerical Recipes. This 
requires an initial guess for the parameters, which we 
find by scanning only the grid points of (j){x, y) (without 
evaluating the interpolating polynomials) for the values 
of the parameters that minimize x^. Then we iteratively 
minimize over the continuous functions, typically arriv- 
ing within less than 10~^ m of the minimum after six 
steps. 



C. Phase Measurements 

We can improve our measurement of the betatron 
phase advance between BPMs by incorporating our 
knowledge of the nonlinear button response. In this mea- 
surement, the beam is excited to small oscillations around 
its equilibrium position (a;o,yo)- Let the phase and am- 
plitude of the AC signal on the «*^ button be represented 
by the complex number C^, and let the phase and am- 
plitude of the horizontal and vertical components of the 
oscillatory beam motion be represented by complex num- 
bers A^; and A^, respectively. To first order, their rela- 



tionship is given by 



(20) 

(21) 
(22) 



The 0i are the functions described in the previous sec- 
tion. Their derivatives are easily calculated from the co- 
efficients of their interpolating polynomials. 

Given the measured Ci, we calculate A^; and A^ by 
minimizing 





Ci = ri^^A^ -f Ti 


yAy, 


where the r^^ 


(x,y) are given by 
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(2:0, an 



4 ^ 



(23) 



Since the ai depend on the closed orbit deviation and the 
values of A.x,y: the minimization must also be performed 
iteratively. The horizontal and vertical phase advance is 
then given by the complex phase of A^ and Ay . When- 
ever a horizontal excitation creates a vertical amplitude, 
or vice versa, this method is used in CESR to compute 
the coupling coefficients also. 



IV. RESULTS 

Testing the new system presents a challenge in that 
we can only produce controlled large amplitude orbits 
with the electrostatic separators. Since the separators are 
calibrated from BPM measurements, they do not provide 
an independent check on our ability to measure large 
amplitudes accurately. Our strategy, therefore, must be 
to use other measurements to check the accuracy at small 
amplitudes, and then confirm the expected linear relation 
between the separator strength and the beam position at 
large amplitudes. 
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FIG. 5; Beam position at various detectors showing little or 
no bunch length dependence. 

To perform a two-dimensional approximation, we ar- 
gued that the bunches are sufficiently long. To verify that 
assumption, we have looked experimentally for a bunch 
length dependence in large amplitude orbits. With the 
pretzel at its nominal value of about 1.5 cm closed orbit 
deviation, the bunch length was calculated from the mea- 
sured synchrotron tune, which we adjust by changing the 
RF accelerating voltage. As Fig. |S1 illustrates, the beam 
position shows little or no dependence over the range of 
bunch lengths we expect in CESR. 

Changing the RF frequency in CESR changes the beam 
energy, and in dispersive regions, changes the beam po- 
sition by up to a few millimeters. Measuring the beam 
position at many different energies allows us to measure 
the dispersion, which we compare to the theoretical value 
from the lattice in Fig. |3 This agreement verifies the 
small amplitude, or linear part of our nonlinear models. 
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FIG. 6: Measured and calculated dispersion. 
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FIG. 7: Beam position at two detectors calculated with the 
nonlinear and linearized methods. 



To observe the large amplitude accuracy of the new 
system, we rely on the electrostatic separators to change 
the orbit amplitude linearly. By increasing the horizon- 
tal separator strength, we observe in Fig. |7| that the or- 
bit calculated with the nonlinear method does show the 
correct behavior, while the orbit calculated with the lin- 
earized formula shows the deviation that was illustrated 
in Fog. 

To demonstrate improvement in two dimensions, the 
voltages on individual horizontal and vertical separators 
were scanned over a regular grid. The measured beam 
positions should also lie on a regular grid, which is shown 
in Fig.|Sl Some sheering is evident in the plot, which may 
be due to coupling of the vertical and horizontal motion 
between the separator and the BPM, or to a rotation of 
the BPM. The pincushion effect is notably reduced with 
the new calculation. 

We use betatron phase measurements to correct the 
difference between the physical optics and the values in 
our model lattice. Without the nonlinear correction, 
large closed orbit distortions hindered this process since 
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FIG. 8: Separator scan. Orbits at detector 9W calculated 
using linearized (dashed) and nonlinear (solid) methods. 
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FIG. 9: Difference in horizontal betatron phase advance be- 
tween data and model with large closed orbit distortion af- 
ter using a phase correction algorithm based on the linear 
(dashed) and nonlinear (solid) BPM evaluation. 



the data we sought to fit did not correspond to the actual 
phase. Figure shows the drastically improved agree- 
ment we can achieve between the model phase and the 
data when the new BPM calibration is used. 



V. BPM CALIBRATION 

The response of a particular BPM may differ from that 
of the computational model (linear or nonlinear) for a va- 
riety of reasons. The leading candidate for this effect is 
the variation in insertion depth of the individual buttons 
(i.e., the distance from the button surface to the surface 
of its cylindrical housing) . This manifests itself as differ- 
ent gains for the signals from different buttons. 

Following the method of 0, Q i we determine the gain 
for each button from the capacitive coupling between 
each pair of buttons. If the ideal coupling between two 



buttons is given by Uij , then the measured coupling will 
be Uij = bibjUij where bi,bj are the effective gains of 
the input and output button, respectively. Symmetric 
pairs of buttons have equal ideal coupling, so U12 = t/34, 
Ui3 = C/24, and C/14 = C/23- 

Using this symmetry for the six measurements Uij {i = 
1, . . . , 3, j = i + 1, . . . ,4) we can calculate the four bi up 
to an multiplicative factor. Normalizing to 61 gives 



61 = 1 
h = \l 



U23U2A 
U13U14 



'U23U 
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(25) 

(26) 

(27) 



These gain coefficients are used to correct the button 
signals before calculating the beam position. 
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FIG. 10: Fractional improvement in the x^ oi the beam po- 
sition fit due to the calibration coefficients ( — 1 = 100% im- 
provement). 

With data drawn from approximately 3700 individual 
beam position measurements, Fig. 1101 shows that when 
these coefficients are employed, the x^ of the fit between 
the measured signals and the modeled signals is signifi- 
cantly reduced. The resulting correction to the calculated 
position is shown in Fig. llll to be approximately 0.5 mm. 



VI. CONCLUSION 

Two-dimensional, electrostatic models of BPM pickup 
response have been used with great success at CESR to 
measure beam position and betatron phase advance for 
large closed orbit distortions. 

Calibration of our BPMs has reduced measurement er- 
rors due to button misalignments. 



1200 



1000 



800 



600 



400 



200 




VII. ACKNOWLEDGMENTS 



FIG. 11: Correction in the calculated position due to the 
calibration coefficients. 



The authors wish to thank REU student Beau 
Meredith for his cahbration of most of CESR's BPMs. 



[1] R. E. Shafer, in Physics of Particle Accelerators, edited by 
M. Month and M. Dienes (American Institute of Physics, 
New York, 1992). [5] 

[2] P. Bagley and G. Rouse, CBN 17 (1991). 

[3] D. Sagan, R. Meller, R. Littauer, and D. Rubin, Physical [6] 
Review Special Topics - Accelerators and Beams 3 (2000). [7] 

[4] S. Krinsky, in Frontiers of Particle Beams; Observa- 



tion, Diagnosis, and Correction, edited by M. Month and 

S. Turner ( Springer- Verlag, Berlin, 1989). 

J. Cuperus, Nuclear Instruments and Methods 145, 219 

(1977). 

G. Lambertson, LSAP Note-5 (1987). 

J. Keil, Ph.D. thesis, Universitat Bonn (2000). 



